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Abstract 


In fusion reactions, the Coulomb barrier selects particles from the high-momentum part of the 
distribution. Therefore, small variations of the high-momentum tail of the velocity distribution can 
produce strong effects on fusion rates. In plasmas several potential mechanisms exist that can pro¬ 
duce deviations from the standard Maxwell-Boltzmann distribution. Quantum broadening of the 
energy-momentum dispersion relation of the plasma quasi-particles modifies the high-momentum 
tail and could explain the fusion-rate enhancement observed in low-energy nuclear reaction exper¬ 
iments. 
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I. INTRODUCTION 


Many-body collisions broaden the relationship between energy and momentnm of qnasi- 
particles: a distribntion of momenta, which can have long tails, characterizes a qnasi-particle 
with a given energy. Therefore, the momentnm distribntion can be very diferent from the 
one obtained using a sharp correspondence between energy and momentum 
typical environments where this effect can be important. 


If. Plasmas are 


Fusion processes select high-momentum particles that are able to penetrate the Coulomb 
barrier and are, therefore, extremely sensitive probes of the distribution tail 12,0,0,1^. 
This broadening of the interacting particle energy-momentum dispersion relation has 


been proposed recently 




as a possible explanation of the strong enhancement of the 


observed low-energy rate of the reaction d{d,p)t in deuterated metal target 

In this paper we study the details of this quantum broadening effect using a simple 
and effective expression for the distributions. In particular, we determine the region of the 
distribution responsible of the effect. Our method is applied to the specihc case of the 
enhancement observed in the d{d,p)t reaction rate. 


II. CHARGED PARTICLE DISTRIBUTION IN PLASMA 


We consider two species (1 and 2) of interacting charged particles with mass, velocity, 
momentum, energy and density: mi,2 , vi,2 , pi,2 , -^1,2 , '^•1,2 • Their fusion reaction rate is 
r = (l + ^ 12 ) ^nin 2 {(TVrei) where the reaction rate per particle is: 


(aVrel) = / d^Pi / d^Pi <hi(pi) <F 2 (P 2 ) aUrei 


( 1 ) 


$ 1,2 (p) are the momentum distributions of particles 1 and 2 and Vrei = |vi — V 2 I is their 
relative velocity. 

The charged-particle fusion cross section a is conveniently expressed as 





( 2 ) 


where S{ep) is the astrophysical factor as function of ^ with p the reduced 

mass and Eq = 2/ic^(ZiZ|Q;7r)^ the Gamow energy. Note that the cross section depends 
on the relative momentum p (to stress this point we write e^); only in special cases there 
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is sharp relation between the energy of the particle E and e^, for instance for free particles 
E = epis the kinetic energy in the center of mass. 

In general we may assume a relation of the form 

1 7 


Q: 


^7 Cp) 


(3) 


ti{E - ep)2 + 72 

where the width 7 = hvcoii depends on the collision frequency Vcoii = nacoiiVcoii- 

Even when the energy distribution is Maxwellian oc exp{E/kbT) (we set the Boltzmann 
constant kb = 1 ), the resulting momentum distribution 


<h(ep)(iep = 


dE d^{E, ep)e 

/n°° P^dp /°° dE 6^{E, ep)e-^/^ 


(4) 


;o ^ JO 

can be non-Maxwellian. We consider the case of a Maxwellian energy distribution, which is 
relevant for the deuteron distribution in metals that is discussed in the next section. A more 
general Fermi distribution, relevant for high-density environments, yields analogous effects 
and, in particular, a power-law tail for the momentum distribution. 

For the sake of concreteness let us consider a Coulombian collisional cross section, acoii = 
e^/e^; the resulting dispersion-relation width is 

2 


2E 


7 = hn-\ — = — 


m 


E, 


X 





(5) 


where the collisional velocity Vcoii = \/2E/m has been used, n is the density and m the mass 
of the colliding particles, deuterons in the present case. For convenience we have dehned the 
energy scale 


Es= — 


f n\ 




m 


2/5 


- Eo = 

no J 


mp\ 1/5 / n\ 
no) 


m 


2/5 


3.02649 eV 


( 6 ) 


where rrie and irtp are the electron and proton masses, Eq = {l/2)a^m(.(? is the Rydberg 
energy, no = (2ao)~^ = 0.843542 x 10^“^ cm“^ a reference density with oq the Bohr radius. 
Then the dispersion relation can be written as 

1 


6 p) 


{ep/Esf^/W^s 


(7) 


EsT^{ep/Es)\E/Es-ep/Esf + E/Es ’ 
and the momentum distribution is conveniently expressed in terms of the scaled variable 
y = 


$ 


((TF|)V 6 ’^ 


T 


X y^^'^dy / dx- 


a/xc" 


X -f 


-{tr 


( 8 ) 
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with the normalization 


AT-i 




dx- 


X + 



( 9 ) 


This distribution should be compared to the Maxwellian one obtained when 6^{E,ep) = 
S{E - ep) 

i^p/T) dep = dz-^y/ze~'' ( 10 ) 

VTT 

with = ep/T. 

Using the scaled variable e^/the distribution depends only on the adimensional 
parameter T /Es- Both the distribution ‘h(ep) and the thermal mean (aVrei) can be obtained 
numerically. We have found an analytical approximation for the important physical regime 


T<^Es 


( 11 ) 


that is sufficiently accurate and allows a better analysis. In fact, since contributions to the 
integral for energies x > 1 are exponentially suppressed, Eq. (HH) implies that the limit 
T/E 5 0 in Eqs. (jHI) and 0 is well-dehned. In this limit the distribution becomes: 


d>r 


(- 


V(TE|)V6 j 


-I der, = I 


V (rifyw’ 0 ) *!> = 




'0 


X + 


where 


No-^ = 


dyy^/^ 


dx 


xe-^ 7ry2r(13/12) 


= 0.51946 


( 12 ) 


(13) 


y^ + x 3(l + v/3) 

Under condition (HU the distribution depends on the temperature T and density n only 
through the single scale parameter ~ which replaces the temperature T 

of the Maxwell-Boltzmann distribution. 

An estimate of corrections to this limiting behavior can be obtained by considering the 
leading corrections to the normalization integral 


iV-‘ (T/Es) = {TIEsf^ + O ({T/Esf^) 

= 0.519457+ 0A37082{T/Esf^^+ 0 {{T/Esf^^) . 

In the limit of small and large Cp the momentum distribution m behaves like 
\im ^o{y)dy = Noy/^y^^‘^dy (l - y/wy^ + 2y^ + 0{y^)) 

6p — 

y/ndy 


lim ^o{y)dy = Nq' 


22/7/2 


3 15 


(14) 

(15) 

(16) 
(17) 


4 











Approximations that retain the leading and next-to-leading asymptotic behaviors of the 
distribution are 

^i{y)dy = Noy/^y^/^dy ^ (18) 

1 -I- 2 y^ 

<S>,u(v)dy = 

The linear combination of the two approximations that maintains the normalization, 
dy{c^i{y) (1 - c)^niiy)) = dy^oiy), reproduces ^o{y) within 2%. 

The most important lesson we learn from this analysis is that the quantum broadening 
due to plasma effects produces a distribution whose typical scale is instead of T: 

these two scales can be very different. For instance, deuterons, m = 2mp, at the density 
n = 4.38 X 10^^ cm“^ have the energy scale Eg ~ 2.03 eV. Therefore, at temperatures 
Ti = 2.44 X 10“^ eV and T 2 = 0.109 eV, which fulhl condition (ITT|) . the scales of the modihed 
distributions are = 0.770245 eV and = 0.988481 eV, respectively. 

This large shift of the particle distribution towards higher energies is demonstrated in 
Fig. Cl where we show *Fo(fp) from Eq. m compared to the Maxwell-Boltzmann distribu¬ 
tion, Eq. m, for the two temperatures Ti and T 2 ; the target deuteron density of Ref. 
has been used as density n of colliding particles: n = 4.38 ■ 10^^ cm“^ . 


III. THE LOW-ENERGY REACTION RATE. 


The d(d,p)t fusion reaction rate has been recently measured using deuterated metal tar¬ 
gets in the 4-20 keV energy range 0, Q, Q, Q| . At low energy these experiments have 
found a considerable higher reaction rate than the corresponding one measured using gas 
targets. Low-energy enhancements are usually explained in terms of electron screening; 
however, the electron screening potential Ue that would reproduce these measurements in 
deuterated metals is of the order of hundreds of eV: this potential is much higher than the 
adiabatic estimate for the maximal screening potential Ue < 28 eV. 

The thermal motion of the target atoms is another mechanism capable of increasing 
the reaction rate; however, the Maxwellian momentum distribution at the experimental 
temperatures gives a negligible effects P, [l^ . 

The observation that large enhancements have been observed in deuterated metals 
but not in insulators |l0l has suggested a possible explanation based on effects of 
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the plasma of electrons in the metal 


HQ 


This simplihed model with quasi-free va¬ 


lence electrons predicts an electron screening distance of the order of the Debye length 
Roeb = \/kbT / {4:7mes{Zey), where Ues is the effective density of valence electrons that can 
be treated as quasi-free. This approach reproduces both the correct size of the screening 


Q,Q. 


potential Ue and its dependence on the temperature: f/g oc 

The problem with this explanation is that the resulting Roeh^ for fhe actual experimental 
conditions, is about ten times smaller than the Bohr radius Oq; the mean number of quasi- 


free particles in the Debye sphere iV^ogft, the so called Debye number is, therefore, much 
smaller than one: Noeh = 'riefr(47r/3)i?pg(j ~ (47r/3)nefT(ao/10)^ ~ 3 ■ 10“^. The picture of 
the Debye screening, which should be a cooperative effect with many participating particles 
(i?Deb should be at least greater than the Wigner Seitz radius, which is of the order of the 
Bohr radius), seems not to be applicable and the observed increase of the d(d,p)t reaction 
rate still missing a consistent explanation. An additional technical inconsistency in the 



Debye screening explanation [10, [lOI is the use of a non-degenerate formula for the screening 
radius in a situation where the electrons are degenerate. 

In this context, we apply the analysis presented in the previous section to the discussion 
of a recent interesting tentative explanation of these puzzling experimental rates, which is 


a0 


based on the quantum-tail effect 

Since it is a good approximation of the experimental situation to consider the pro¬ 
jectile distribution monoenergetic with a sharp relation between energy and momentum 
S^{Ebeam,^p) = S{Ebeam “ Cp), the relative velocity is function only of and of the angle 
■d between pt and the beam: = 2{Ebeam ep^ — 2 ^ Ebeam^pt cos -d) / m. Following the 

analysis of the previous section, target particles have a Boltzmann-Gibbs energy distribu¬ 
tion with a relation between energy and momentum that is broadened by plasma quantum 
effects and it is given by Eq. ©. The effective momentum distribution of the target par¬ 
ticles is, therefore, not Maxwellian but given by Eq. (jSl). Substituting this distribution of 
target particles, the sharp monoenergetic distribution for projectile particles, and the above 
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relation for the relative velocity in Eq. (CD the reaction rate per particle becomes 


{crVrel) 


^ '^rel 1 


^ ^rel 


^ + 1 POC poo 

27rm^^ / dcosi} / dE / (iepj-y2e^5^(E, e^J 


dcos^ / depM^pJ{TElY/^,T/Es)aVrei , 


'-1 


TTId 'i^rel 





(20) 

( 21 ) 


where Ecm = and p = m/2. Since condition T -C is verihed, we can use the 

simplihed form in Eq. m 

p1 poo 

{(^Vrei) = ^ j dcosd j dcp ,$0 (cp,/(TE|)^/®)cr ■ (22) 

As we have analyzed in the previous section, the momentum distribution resulting from 
the quantum broadening with the chosen collisional cross section has two main features: a 
peak at energies of the order of instead of T and a power-law tail that decreases 

as Cp ' instead of the exponential cut-off. It is, therefore, physically interesting to separate 
the contributions from the peak, the high-momentum tail and the low-momentum part of 
the distribution, so that we can understand which feature(s) of the modihed distribution 
give(s) important corrections to the rate. To this purpose we dehne the peak of the target 
momentum distribution as 0.54(TE'|)^/® = ei < ep < = 1.15(TE|)^/®: this region includes 

about 50% of the particles and we call the corresponding contribution to the reaction rate 
per particle {avrei)cl the contributions from the low- (cp < ep about 10% of the particles) 
and high-momentum (e/j < Cpi about 40% of the particles) parts are indicated with the 
subscript L and H. 


{aVrel) = {crVrel)L + ^rel) c + {(^ V^el) 


/H 


(23) 


We have used S{Ecm) — 5'o = 43 keV b (the error is < 6% for Ebeam < 10 keV): then 
the d integral can be done analytically in terms of the incomplete Gamma-Euler function; 
the remaining integrations have been performed using the Gauss-adaptive method. 

The resulting rate is shown in Fig. |21 as the thin solid line. In the same hgure are 


shown for comparison the experimental data 


and the other theoretical curves: the thick 
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solid line shows the rate for bare nuclei, while the dot-dashed line shows the rate with 
electron screening in the adiabatic limit, which should provide an upper limit to the screening 
potential, /g = exp j. 

As it is apparent from Fig. |2l the quantum-tail effect is in fact capable to produce a 
strong enhancement of the reaction rate, but this effect starts only below Ebeam ~ 2 keV; 
on the contrary the experimental excess starts at energies two or three times higher. 

From the analysis of the separate contributions of the three regions (low, central or 
peak, and high) of the target momentum distribution (see Fig. El), we observe that the 
increase of (avrei) at low energies {Ebeam ~ 2 keV), shown in Fig. |21 is caused mainly by 
the high momentum particles in the power-law tail, the {(TVrei)fj term: the peak and the 
low-momentum region seem not contribute to this increment in the present situation. 


IV. CONCLUSION 

We have studied the effects of the quantum broadening of the relation between energy and 
momentum due to a specihc collisional cross section. An important effect is the considerable 
modihcation of the resulting momentum distribution: 

(1) the central part of the distribution is shifted from T to (TEgY^^ where the energy scale 
Es, Eq. ©, grows as with the density n; 

(2) the high-momentum tail decreases as a power Cp ' instead of having an exponential 
cut-off (see Fig. [T]). 

We have applied this quantum-tail effect to nuclear fusion processes between charged 
particles at sub-barrier energies of the order of few keV and compared our results with the 
experimental data relative to the d(d,p)t reaction with deuterated target. 

Our calculation shows that the quantum-tail effect produces an important increment 
of the observed reaction rate enhancement at very low energies (~ 2 keV). However, this 
mechanism cannot reproduce the experimental rate for deuterium, which has been found to 
increase already at higher energies (~ 6 — 8 keV), as shown in Fig. |21 

We have also analyzed more in details the effect of the modihed momentum distribu¬ 
tion on the reaction rate by breaking up the contributions from target particle in three 
regions: the low-momentum, the central or peak, and the high-momentum region. The 
strong enhancement of the rate is due essentially to the particles in this last region: the 
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high-momentum power-law tail of the distribution. 

We are extending our results by using other collisional cross sections and investigating 
the temperature dependence of the mechanism. 
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cVei (cm^sec"^) S^^p(keVb) 




FIG. 2: Reaction rate per particle {(JVrei)T=o panel) and the corresponding astrophysical factor 
S = {crVrei)T=o (^\/(right panel) for the d{d,p)t reaction as function of 

the beam energy. The experimental data ^ are compared with three theoretical curves: bare 
nuclei (dashed line), screened nuclei with the adiabatic potential 14 = 28 eV (dot-dashed line), 
and our calculation that includes the quantum-tail thermal effect (solid line). 
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FIG. 3; Contributions to the reaction rate per particle (ct Vrei) coming from the three different regions of 
the target momentum distribution: low-momentum region (dash-dotted line), central or peak region (dashed 
line), and high-momentum or tail region (solid line). 
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